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Abstract 

The complete set of 0{a) QED corrections with soft photon exponentiation to the process ep —>■ 
eX in mixed variables {y = yh, = Qf) is calculated in the quark parton model, including 
the lepton-quark interference and the quarkonic corrections which were unknown so far. The 
interference corrections amount to few percent or less and become negligible at small x. The 
leading logarithmic terms proportional to InlQ"^ /m'^) from radiation off quarks are discussed and 
the non-logarithmic quarkonic corrections found to be negligible for almost all experimentally 
accessible x and y. 
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1 Introduction 



Traditionally, deep inelastic ep scattering, 

e{h)+p{p^)^e{h)+X{p2), (1) 

accompanied by radiative corrections from the reaction 

e(fci) + p{pi) e{k2) + X{p2) + n^{k), (2) 

was studied in terms of kinematic variables, which could be determined exclusively from the 
electron. The new detector generation at HERA gives additional access to the hadron kinematics. 
For the radiative process (H) the 'leptonic' variables differ from those, which are determined 
partly or exclusively with information about the hadronic kinematics. QED corrections from the 
diagrams of figure |l| (henceforth called leptonic corrections) with soft photon exponentiation have 
been determined in a model independent approach in several variables [|I], 0. The same has been 
undertaken in leading logarithmic approximation, including higher order corrections p|, |^, for a 
larger variety of variable sets [Q. For a complete 0{a) calculation, one has to use the quark- 
parton model and include corrections from the diagrams of figure ^ and from their interferences 
with those of figure 0. In terms of leptonic variables, this has been done some time ago 

In this letter, we report on the corresponding corrections in terms of mixed variables. These 
variables are defined as follows: 

ym = yh = — 7^ — , Q^ = Qi = [h - ^2) , Xm = — t;- (3) 

In the next section, we explain the structure of the corrections and give explicit expressions 
for them. Section 3 contains a discussion of photonic bremsstrahlung from quarks and section 4 
numerical results. The technically involved calculations will be described elsewhere. 



2 The bremsstrahlung corrections 

The cross section of reaction (2) may be written as follows^: 
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(4) 



(5) 



andy±(y) = l±(l-y)2._ 

The sum over Q and Q extends over the quark content of the proton. The second sum with 
a = e,q,i extends over the photonic corrections originating from the lepton legs, from the quark 
legs, and from their interference. Correspondingly, Ce = Ql,Cg = Q 



QsQq- We use the 

convention = —1. The couplings Qf,Vf,af are collected in overall factors of the vector type, 
V{B,p), and axial vector type, A{B,p) and the running coupling agEj^iQ"^), the Fermi constant 
Gp, and the ratio of the photon and Z propagators in effective coupling strengths K{B,p): 

^If not stated differently, the variables x, y, are mixed variables. 
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V{a,h) = K{a,b;B,pePQ)V{B,p,pQ), (6) 

B=-y,I,Z 

K{a,b;B,p) = XB^ia) XBiib) K{B,p), 

Q2 

^^^^^ = aQ^ + Ml" 

Since the interaction of leptons with quarks proceeds via photon or Z boson exchange, we have 
three types of cross section pieces labeled with B = j, Z, I, where the latter denotes again the 
interference. The [-B1-B2] = 11, Z'j, ZZ correspond to i? = '-y^I^Z. The V{B,p) and K{B,p) 
contain weak one loop corrections; they are defined in The axial couplings A{a, b) are related 
to the couplings A{B,p) of |^ in the same manner as is explained above for the vector couplings. 
The sign factor is p = PePq = +1 for PP, A A scattering and p = —1 for PA scattering, where we 
denote particles with P and anti-particles with A. 

^ Deviating from we took the Xz{^) out of the definition of the K{B,p). An update of the weak corrections 
is described in [^. 
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The functions S^, Sq contain factorizing soft photon corrections and the corresponding vertex 
corrections^: 



Seix,y\ml 



'in ^ - 1 Vn 



1 — a; 



X 



+ 2Li2(l) -Lialxy) + Li2 



3 

+ o 1^ ^ ~ 2 - ln(l - y) ln{xy) 
(1 - x)y 



Sq{x,y\ml) = Se{x,y\ml) 



Li2(l -y) + Li2 



1 — X 



(7) 



X 



With the exception of the last two terms in (||) the expressions for Se{x,y\ml) and Sg{x,y\mQ) 
agreeQ. 

The soft photon interference term Si has to be combined with the 77 and 7Z boxes in order 
to get an infrared finite answer; the latter are added as an extra piece after the divergences are 
canceled, and it is: 



S,ix,y) = ln(l-|/)ln^^-^4^^-6Li2(y) + 2Li2(xi/)-2Li2 



(1 - x)y 
1-2/ 



(9) 



The corrections B^'^ and B^'"^ in are due to real photonic bremsstrahlung. Although we 
leave out all details of the involved calculations, a remark on the phase space parameterization 
should be made. The double-differential bremsstrahlung corrections are the result of a threefold 
integration, which has been organized such that two of the integrals may be performed analytically. 
At the parton level, there are two nontrivial integrals to be performed whose result has to be 
convoluted with the parton momentum distribution: da^^'^"^ = a^/{27iS) ■ J2q,q fQ{xh){dxh/xh) 
dQ'jdQ\dyidT k\M^^q^\'^ ■ We chose the following parameterization of the phase space: 
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(10) 



with Q\ = {p2 — PiY- The integration variable z may be larger or smaller than 1. As may be 
seen in figure 39 of for these two cases the physical regions for the integration over yi (at fixed 
values of Q\ or, equivalently, of z) are different, which in turn leads to different integrands in (^ 
for the subsequent integration over z. The integral over yi and that over the photon phase space, 
indicated by dTk, are performed analytically. In fact the latter one is transformed into an integral 
over dz2 / {^y/Wz) , where Z2 = —2kk2 is the denominator of the final state electron propagator 
before an emission of a photon . 

The hard bremsstrahlung functions in (H) are: 



B„ ( 
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(12) 



The exponentiation of the soft photon term (InQ^/m^ — l)ln((l — j/)(l — x)/x) is performed only for the 
leptonic correction 5*6 Q . 

^ Here we treat the quark masses as independent input parameters. Instead, one could use the relation ruq = 
XhMp, where Mp is the proton mass as it is suggested by kinematics. Then, the square bracket in (||) would change 
the sign. In integrated expressions like (p|), it is then — XmMp. 
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The bremsstrahlung vector couplings are: 

K = c,V{z,z), V, = c,V{l,z), V, = c,V{l,l). 



(13) 



The bremsstrahlung axial couplings are defined analogously. The different arguments of V{a, h) 
refiect the different momentum flows in the exchanged boson propagators in flgure |1| and 
figure I (Q2). 

For a = e,i,q, the n runs from 1 to 2,1,2. The functions Fan are multiplied by a factor 
1/(1 — 2;) which becomes singular at 2; = 1 due to the infrared singularity of the corrections in this 
phase space region. This singularity is regularized by a subtraction: 
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The leptonic corrections differ by a simple factor: 
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The lepton quark interference corrections have no mass singularities: 
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The integrands in the two integration regions in (HI) are related by a symmetry relation: 
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X — Le, Lj, Wj, Fq, Lq, Uq. 



(22) 
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Finally, we have to discuss the last term in (^), d?a'^°^ / dxdy. This is the finite part of the 
contribution from the interference of the 77 and 7Z box diagrams with the Born diagrams. It 
contributes to the factorizing part of the lepton quark interference. Its variables are defined from 
the Born kinematics. For this reason, we may refer to eqs. (2.1), (3.1), (3.4) and (C.ll), 
(C.14). The latter two formulae contain the corrections from the box diagrams: Si{B; a,b)''°^ = 
Si{B; a, h) - S'°^\ B = -f,Z, where it is = Aa"^ \n{b/a) ln[(l - x)/x]. 



3 Photonic bremsstrahlung from quarks 

In this section, we discuss the connection of our complete 0{a) results with expressions to be 
expected from a LLA treatment. For this purpose, we rewrite the terms proportional to ln(Q^/mQ) 
in the first of the integrals in (§) and perform a change of integration variable, z —>■ 1/z: 
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In the second integral, one has to collect terms correspondingly (here without change of integration 
variable): 
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Further, the soft photon correction (^) has to be rewritten: 
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As is indicated by the choice of indices, the selected corrections may be described by the 
following generic expression for the collinear radiation of photons from the initial or final state 
quark: 
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where a = ini, fin and niQ is a constant quark mass. The appropriate choices of initial and final 
state scalings are 



im 

Pi = ZPl, 



fin 1 
P2 = -P2, 

Z 



and the corresponding Jacobian is 



d{x,y) 



d{x,y) 



(35) 



(36) 



Further details about scaling of kinematic variables may be found in tables |^ and These tables 
contain, in addition, the scaling properties for other experimentally interesting sets of kinematical 
variables. 

We now may address the following question: Is it possible, as has been proposed for the case of 
leptonic variables 0, ||, , to absorb the quark mass dependences completely by a redefinition of 
the quark distribution functions? For initial state radiation, this is evidently the case. The LLA 
cross section in mixed variables has the same structure as that in leptonic variables and may 
be treated in the same way. Otherwise, the contribution would be substantially overestimated (see 
[^). This holds also for the other variables in table 1 with the notable exception of the hadronic 
ones. 

In leptonic variables, there is in accordance with the KLN theorem no LLA contribution from 
final state radiation |^ . For mixed variables (and other ones which are also not totally inclusive 
in the hadronic final state (see table 2)), this is different: There is a LLA contribution and its 
Born function -Bo(z/, z) in ( ^9]) does not factorize, preventing the framework as developed in Q to 
be operative here. 
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Table 1: The definition of scaling variables for the leading logarithmic corrections due to photon 
emission from the initial quark leg. Those variable sets of table 2, which are not contained here 
transform as the leptonic variables do. The definitions of variables are explained in detail in [3,4]. 
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Final state radiation 
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Table 2: The definition of scaling variables for the leading logarithmic corrections due to photon 

emission from the final state quark leg. Variable sets without rescaling have no final state LLA 
radiation (leptonic and double angle kinematics). The definitions of variables are explained in 
detail m [3,4]. (*) zo = [l- 2xy{l -y)-{l- Axy{l - y))^/^]/[2x{l - yf] 
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4 Numerical results 



The QED corrections will be described by the following correction factor: 

^mi^ = d2^Born X 100% ^ 5, + S, + S,. (37) 

QED corrections in mixed variables are quite different from those in leptonic variables. Numer- 
ical results have been obtained with HECTOR 0, a Fortran program for the calculation of QED, 
electroweak, and QCD corrections to deep inelastic lepton nucleon scattering in different vari- 
ables and different formal approache^: LLA - leading logarithmic approximation, with possible 
inclusion of higher order corrections; MI - complete 0{a), plus (or without) higher order LLA, 
[quark-parton] model- independent; QPM - complete 0{a), plus (or without) higher order LLA, in 
the quark-parton model. The latter approach is advocated here. 

As a matter of fact we mention that the leptonic corrections 6e in mixed variables determined 
here in the QPM agree numerically exactly with those determined in the MI (when using the same 
structure functions). The latter rely on completely different expressions and are described in 
where also the general features of QED corrections in mixed variables are explained. 

On top of the leptonic corrections 6^, which are the numerically largest ones, the effect of 
the interference corrections 6i is shown in figure |^. They are positive at intermediate values 
of y and become negative at y > 0.7, nowhere exceeding in magnitude 1% except very large 
values of x and y. The additional influence of the quarkonic corrections is shown in figures | 
and ^ In figure ^ we present the effect of the quarkonic corrections, calculated for definiteness 
with the assumption mq = XhMp. Their effect is quite sizeable due to the unphysical quark mass 
singularities. Excluding both the initial and final state mass singularities from the prediction (and 
assuming now the case of constant quark masses, which chooses the sign discussed in footnote 4), 
we get figure |^. In this case the quarkonic corrections are rather small and the situation as shown 
in figure ^ is re-established. 

The exclusion of initial state mass singularities is proposed in 0. There it is shown that the 
photonic initial state LLA corrections, although modifying the QCD evolution equations, give 
rise to negligible numerical effects in the kinematical domain of HERA. The exclusion of final 
state mass singularities might be motivated by experimental reasons^. It is scarcely possible 
to distinguish between bremsstrahlung photons and those from 7r° decay, say. For an adequate 
treatment one has to redefine the kinematical variables by adding the energy of a final state 
photon to the hadronic energy. In such variables, the final state LLA corrections would vanish in 
accordance with KLN theorem. 

To summarize, the QED corrections in mixed variables have a substantially different behavior 
compared to those measured in leptonic variables. The bulk of the corrections comes from the 
leptons. If the data will reach an accuracy of the order of a percent one has to take into account 
also the interference bremsstrahlung corrections which, together with the quarkonic ones, have 
been calculated here for the first time. The major part of the latter is associated with quark mass 
singularities and may be excluded from the corrections. 
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Figure 3: Comparison of the radiative corrections S^/^ — 5^ + 5i (broken line) with the leptonic 
corrections Sg (solid line) at HERA in mixed variables. 
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Figure 4: Comparison ofS^^l — 5e + Si (broken line) with Smix — Se + Si + 5q (solid line) at HERA. 
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